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The singularity structure of charged spherical collapse is studied by considering the evolution of the gravity- 
scalar field system. A detailed examination of the geometry at late times strongly suggests the validity of the 
mass-inflation scenario [B|. Although the area of the two-spheres remains finite at the Cauchy horizon, its 
generators are eventually focused to zero radius. Thus the null, mass-inflation singularity generally precedes 
a crushing r = singularity deep inside the black hole core. This central singularity is spacelike. 
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PACS: 04.40. Nr, 02.60.Cb, 97.60IT, 04.20Dw 

The radiative tail of gravitational collapse decays as an 
inverse power of time leaving behind a black hole char- 
acterized only by its mass, charge and angular momen- 
tum [[[] . The simplicity of the external field at late times 
contrasts with the highly dynamical interior which lies 
beneath the event horizon 

It is generally accepted that the "tunnel" through a 
charged or rotating black hole is destroyed by pertur- 
bations which propagate into the hole after its forma- 
tion [EjU^]. The manner in which this occurs has been 
studied in great detail using simplified models ||.@,|. 
and the results suggest that the Cauchy horizon (CH) in- 
side a charged (or rotating) black hole is transformed into 
a singularity at which the "Coulomb" component of the 
Weyl curvature diverges. The null generators of the CH 
contract slowly under transverse irradiation, so it is an- 
ticipated that this null singularity eventually gives way 
to a central spacelike singularity deep inside the black 
hole core. It should be noted that this spacelike singu- 
larity is not expected to precede the CH as some authors 
have argued Indeed, geodesic observers falling into 

the black hole at late times will generally encounter only 
the null singularity. 

Some disquiet exists regarding this picture of the 
black hole interior. Yurtsever [|| has argued, on general 
grounds, that complete destruction of the CH should be 
expected once generic perturbations are considered. He 
bases his discussion on experience with colliding plane 
wave spacetimes — he has shown that plane wave Cauchy 
horizons are replaced by spacelike singularities in the 
presence of generic plane-symmetric perturbations [ p~0| | . 
Gnedin and Gnedin 11, on the other hand, have per- 
formed a numerical integration of the spherical Einstein- 
Maxwell-scalar field equations. They established the ex- 
istence of a spacelike singularity inside a charged black 
hole coupled with scalar matter. Their analysis stressed 
the behaviour of the central (r = 0) singularity, however 
they did not consider the possibility that it intersects the 
Cauchy horizon. 

This letter reports on an independent numerical in- 
vestigation. We find evidence that the analytic models 



demonstrate the essential features of black hole internal 
structure. 

The general spherical line element can be written as 

ds 2 = -ggdv 2 - 2gdvdr + r 2 (d8 2 + sin 2 Odcf?) , (1) 

where g = g(r,v), g = g(r,v) and v is advanced time. 
We solve the field equations G a p = 8ir(E a p + T a fj) where 



£% = (e 2 /87rr 4 )diag(-l, -1,1,1) 



(2) 



is the standard electromagnetic contribution to the 
stress-energy, and 



>,a<P,f3 



TjQaP (</>,7</>' 7 ) 



(3) 



is the stress energy for a massless, minimally coupled 
scalar field (f>. It is convenient to write h = <f> and then 
introduce the scalar h(r, v) by 



(rh), r = h . 

The field equations can now be written as 

(lng) tr = r- 1 (h-hf , 
(rg), r = 5 (1 - e 2 /r 2 ) , 

while the wave equation 0(f> = becomes 
h. v 



^h r = ±(h-h)(g[l-e 2 /r 2 ]-g) 



(4) 

(5) 
(6) 

(7) 



The pioneering numerical integration of these equa- 
tions, with e = 0, was performed by Goldwirth and Pi- 
ran jyj. We have employed a similar algorithm. Us- 
ing the method of characteristics Eq. (0) is recast as 
a set of 2n ordinary differential equations, where n is 
the number of radial grid points. The radial integrals in 
Eq. (^)-Eq. (||) are discretized according to the trapezium 
rule, while we use a Runge-Kutta scheme for the ordi- 
nary differential equations. The resulting code is locally 
second order accurate, and has been tested on Reissner- 
Nordstrom spacetime and the exact self-similar solutions 
in |12|. The details will appear elsewhere. 
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To aid the comparison to previous work |3|,|4| we intro- 
duce the mass function 

m(x a )= r -(l + e 2 /r 2 -g/g) . (8) 

Moreover the sole non-vanishing Newman-Penrose Weyl 
scalar 

-tf 2 = \c H h = [m(x a ) - e 2 /r}r- 3 , (9) 

is simply expressed in terms of this mass function and 
provides a direct measure of curvature near the CH. 

Characteristic initial data are supplied on a pair of 
intersecting null hypersurfaces as shown in Fig. 1. The 
spacetime is Reissner-Nordstrom when v < vq. After 
this advanced time we consider the evolution of a charged 
black hole with infalling scalar matter. Since our primary 
interest is the black hole interior it is useful to choose the 
outgoing null hypersurface, T, coincident with the event 
horizon. This is achieved by first setting g\r = — 1, then 
specifying r\r such that 

lim r|r = constant . (10) 

v — >oo 



r= 




® 



FIG. 1. A spacetime diagram showing the setting of the 
numerical integration. The spacetime is Reissner-Nordstrom 
for v < vq. Scalar field falls into the black hole across the 
event horizon, F. PEH is the past event horizon, located at 
i) — — oo. AH is the outer apparent horizon, and I AH is the 
inner apparent horizon of the charged black hole. A couple 
of lines of constant r are shown in light gray. The Cauchy 
horizon (CH) is a singular hypersurface which contracts to 
meet the central singularity at r = 0. All singularities are 
indicated by thick lines in the diagram. 

The dynamical equations are supplemented by 

one further equation along each outgoing null hypersur- 
face: 



dg _(d\ng 1 r 2 2 _A ( d(j>\ 

Tv= g \—-^ 9{r - e) - g ^) + 2r \Tv) ■ 

(11) 

This is Raychaudhuri's equation for the outgoing null 
rays, and is conserved by the evolution. Using the char- 
acteristic equation dr/dv = —g/2 to obtain ~g\r, it is 
straightforward to integrate Eq. (|ll]) along T, with the 
appropriate initial conditions at v ~ vq, for </f>|r- 

Beneath the event horizon of the black hole outgoing 
null rays contract rapidly (see Fig. 2). Consequently an 
initially uniform radial grid becomes highly non-uniform 
at late times, leading to significant errors in ~g. We have 
employed a non-uniform grid along the initial ingoing 
null ray to alleviate this problem, making the radial grid 
points more dense in the vicinity of the event horizon. 
A significant increase in (accurate) integration time is 
obtained by this method, and the internal mass function 
[defined in Eq. (||)] reaches values up to about 10 50 times 
the external black hole mass before errors become signif- 
icant. 

We have integrated the equations for several choices of 
initial data along T, and have obtained the same qualita- 
tive behaviour. Table || summarizes the results obtained 
for the two initial data sets characterized by 

, f r + - [3v~ v , . 

r r = { a i 2^ (12) 

' [ r+ — p exp (— pv z ) v ' 

where (3 > and p > are real parameters. It is conve- 
nient to fix dimensions such that the asymptotic Bondi 
mass is unity, thus the radius of the event horizon ap- 
proaches r + = 1+vl — e 2 as v — * oo. The inverse power- 
law data along T is taken as representative of the physi- 
cal situation where the flux of radiation across the event 
horizon exhibits such a decay Jlj| . We have included the 
second data set to emphasize that the non-linear instabil- 
ity of the inner horizon is present even for perturbations 
having compact support on the event horizon. This re- 
sult should be expected based on linear analysis of the 
problem 0. 

TABLE I. Late time behaviour of the metric along char- 
acteristics labelled by u. Multiplicative factors which depend 
only on it have been omitted. The constants 7, /1, and a are 
discussed in the text. 



Initial data 


g(u,v) 


g(u,v) 


m(u, v) 


r+ — f3v~ p 


e -iv 


v -( p +1 + <t) 


u -(p+l + CT) e 7" 


r+ - /3e ( - pl,2) 




e -2„v 





In Reissner-Nordstrom spacetime all outgoing null 
geodesies which originate inside the black hole intersect 
the CH within a finite affine parameter. However, the 
presence of the scalar field modifies the causal structure 
of the spacetime so that some outgoing null rays termi- 
nate at r = in a finite coordinate time v (i.e. before 
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intersecting the CH). This important result is demon- 
strated in Fig. 2, which shows the radius as a function of 
external advanced time along a selection of outgoing null 
rays. In this coordinate system the lightcones tip over at 
r = indicating that the singularity at this hypersurfacc 
is spacelike. 
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FIG. 2. The null rays in the rw-plane. The inner-most rays 
terminate at r — in a finite advanced time, however there 
exist many geodesies which approach a finite radius at large 
v. Subsequent figures show ~g and g along the five outermost 
rays indicated by different line types. The initial data is ex- 
ponential with j3 = 0.26, p = 1 and e 2 = 0.4. 

Generally we have found the outgoing null rays to be- 
long to two distinct classes, those which reach r = and 
those which do not. The latter geodesies are of great- 
est interest to us here and below. Each member of this 
class approaches a fixed radius at late times, indicating 
the existence of a CH in these solutions. Notice that 
the final radius is different for each geodesic; the CH is 
not a stationary null hypersurface, as it is in Reissner- 
Nordstrom spacetime, rather it contracts slowly under 
transverse irradiation. Moreover, the CH is the locus of 
a null, precursory singularity which precedes the crushing 
singularity at the origin (see Fig. 1). Gncdin and Gnedin 
examined only the central singularity in j^], overlooking 
the existence of the null portion. 

One might worry that the cumulative effect of infalling 
scalar field could eventually focus all outgoing null rays to 
zero radius. To allay such doubts we have examined the 
expansion rate ~g = —2dr/dv along outgoing null rays. 
Typically g reaches a maximum along each geodesic, then 
it exhibits a period of exponential decay followed by a less 
pronounced, but definite, approach to zero which is char- 
acteristic of the initial data (see Fig. 3). Table |] gives the 
form of ~g for large v, demonstrating that the assumptions 
used in asymptotic calculations are valid . 

Our initial conditions along T reflect the relaxation 
of the external gravitational field to Reissner-Nordstrom 
spacetime at late times. It seems reasonable that this 



property should continue to hold inside the event hori- 
zon as far as some imprecisely defined radius, where 
non-linear gravitational effects become important. This 
expectation is indeed realized in our numerical simula- 
tions, justifying the assumption used in |D| to analyti- 
cally estimate the effects of the scalar field on the charged 
black hole interior. We find a thick (in terms of r) layer 
where the presence of scalar matter produces only slight 
deviations from an exact Reissner-Nordstrom solution. 
This effect is most pronounced in the function g at late 
times; near to the event horizon g ~ — 1 as in Reissner- 
Nordstrom spacetime, however it rapidly approaches zero 
beyond this region. 
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FIG. 3. These plots show In \g\ along a selection of outgoing 
null rays which intersect the CH. The top graph shows In \g\ 
against In \v\ for the inverse power data in Table I. The late 
time fall-off is clearly a power law too. The graph is In \g\ 
against v for the exponential data. The linear relation in this 
graph indicates an exponential fall-off at late times. The cusp 
in these figures corresponds to a change in the sign of g where 
the outgoing null ray intersects the outer apparent horizon of 
the black hole. 

Figure 4 shows hx\g\ plotted against advanced time 
along selected null rays. At late times g ~ e -7 ", where 
7 depends only on the charge e. We find 7 ~ k_ as pre- 
dicted by analytic models B, where k_ = \/l — e 2 / (1 — 
yl — e 2 ) 2 is the surface gravity of the inner horizon of a 
static black hole with equivalent mass and charge. This 
result holds to within about 10%. 

While a suitable coordinate transformation can render 
the metric non-singular [0], the exponential decrease in 
g is reflected in the growth of curvature as the CH is ap- 
proached. In fact C a ^ 5 C al3l s ^ 48m 2 (u, v)/r 6 {u, v) ~ 
g 2 (u,v)e 2jv along outgoing null rays, (u labels the out- 
going null hypersurfaces on which we have examined the 
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functions.) For a black hole of about one solar mass 
we have been able to follow the evolutions to curvatures 
around Planck levels. 
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FIG. 4. In \g\ against advanced time along the outgoing 
null rays. The asymptotic form g ~ e~ JV is in remarkable 
agreement with predictions based on simplified models. 



we really get? Without the entire solution (up to and in- 
cluding the CH) it is difficult to quantify this. However, 
in terms of a Kruskal-like coordinate V = —e~ K - v which 
goes to zero on the CH we have reached values as small 
as V = — 10~ 150 during our integrations. It might seem 
tempting to choose such an advanced time coordinate ab 
initio, so that the CH is located at a finite coordinate 
distance Vch- However a high price is paid for such a 
choice; in the new coordinates the initial data g\r and 
g\r become badly behaved as V — ► Vch leading to sig- 
nificantly decreased integration times. 

In conclusion, we believe that the effects we have de- 
scribed above are representative of the asymptotic struc- 
ture of the true black hole interior and that a null, mass 
inflation singularity is present along the CH. Further- 
more, the null CH singularity is a precursor of the final 
spacelike singularity deep within the black-hole core. A 
detailed account of this work is in preparation. 

We are grateful to Ian Moss and Werner Israel for use- 
ful discussions. This research was supported in part by 
the EPSRC. 



Analytic models have relied heavily on the results of 
linear perturbation theory to provide information about 
energy fluxes along the CH. Therefore, we have exam- 
ined the flux of scalar field parallel to the CH in our 
solutions, finding qualitative agreement with perturba- 
tive analyses. It is straightforward to see this remarkable 
feature directly from Raychaudhuri's equation (|ll]) and 
the results in Table |. If l a = dx a /dv is a lightlike gener- 
ator of the outgoing characteristics then the flux of scalar 
matter across these surfaces is 



T := AnT aP l a l p = (d<p/dvf 



(13) 



Since r — > constant and \g\ decays exponentially at late 
times, Eq. (O) implies that 



(14) 



T ~ Ag + B dg/dv as v 



where A and B generally depend on the outgoing null 
surface considered. Now, if the flux across the event 
horizon at late times is J-\r ~ v~ p ~ 1 then T ~ T\y v~ a 
with a > 0, along null rays which intersect the CH. The 
value of o~ (lying in the range 0.3 < a < 0.75 for the 
cases we examined) depends on the charge, and differs 
from predicted values 0,|l5| by more than the numer- 
ical uncertainty. Generally, if T\r decays faster than 
exp[— 2k+v] for large v, then T ~ exp[— 2/zu] where fi ~ 
k+ = Vl - e 2 /(l + Vl -e 2 ) 2 (to within less than 10%). 
The striking agreement with predictions of linearized the- 
ory |t]] suggests that perturbative arguments ||[l5| are 
actually valid deep inside the black-hole core. 

The most significant shortfall in our analysis is the lo- 
cation of the CH at v = oo. So, how close to the CH do 
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